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«N : NON-LTE MODEL ATOM CONSTRUCTION 

O 

^ ' Norbert Przvbilla^ 

Abstract. Model atoms are an integral part in the solution of non- 
. LTE problems. They comprise the atomic input data that are used 

Mm , to specify the statistical equilibrium equations and the opacities and 

■ emissivities of radiative transfer. A realistic implementation of the 

(-H ' structure and the processes governing the quantum-mechanical system 

^ of an atom is decisive for the successful modelling of observed spectra. 

We provide guidelines and suggestions for the construction of robust 
and comprehensive model atoms as required in non-LTE line-formation 
computations for stellar atmospheres. Emphasis is given on the use of 
^ ' standard stars for testing model atoms under a wide range of plasma 

conditions. 
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1 Introduction 



> 

o 

Astrophysical plasmas like stellar atmospheres, gaseous nebulae or accretion disks 
' are not in any sense closed systems, as they emit photons into interstellar space. 
' Therefore, the thermodynamic state of such plasmas is in general not described 
I well by the equilibrium relations of statistical mechanics and thermodynamics for 
\ local values of temperature and density, i.e. by local thermodynamic equilibrium 
^-H . (LTE). The presence of an intense radiation field, which in character is very dif- 
^ \ ferent from the equilibrium Planck distribution, results in deviations from LTE 
' (non-LTE) because of strong interactions between photons and particles. The 
r> ' thermodynamic state is then determined by the principle of statistical equilib- 
rium. All microscopic processes that produce transitions from one atomic state 
■ ■ ■ ' to another need to be considered in detail via the rate equations. A fundamental 
complication is that the distribution of the particles over all available energy states 
- the level populations or occupation numbers ~ in turn affect the radiation field 
via the effects of absorptivity and emissivity on the radiation transport. What 
is required is a self-consistent simultaneous solution of the radiative transfer and 
statistical equilibrium equations. 
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A model atom is a collections of atomic input data required for the numerical 
solution of a given non-LTE problem. It is a mathematical-physical approxima- 
tion to the quantum-mechanical system of a real atom, and its interaction with 
radiation and with other particles in a plasma. A model atom comprises, on one 
hand, data to specify the structure of the atom/ion like energy levels, statistical 
weights and ionization potentials. On the other hand, the transitions among the 
individual states need to be described, requiring oscillator strengths, cross-sections 
for photoionization and coUisional excitation/ionization, etc. The number of levels 
in a model atom amounts typically to several tens to several hundred in modern 
work, and the number of transitions from hundreds to many (ten-)thousands. 

As only a very limited amount of atomic data have been determined experi- 
mentally up to now - mostly energy levels, wavelengths and oscillator strengths -, 
most of the data have to be provided by theory. Large collaborative efforts have 
been made to compute the data required in astrophysical applications via ab-initio 
methods. The Opacity Project (OP; Seaton[T987j Seaton et al. ll994p and its suc- 
cessor the IRON Project (IP; Hummer et al. I1993P provided enormous databases 
of transition probabilities and cross-sections for photoionization and excitation via 
electron impact. Many smaller groups and individuals have contributed additional 
data, most notably Kurucz (see e.g. Kurucz I2006P in a tremendous effort lasting 
already for about three decades. Ab-initio data for radiative processes between 
levels of principal quantum number n < 10 are available for most of the ions of 
the lighter elements up to calcium, and for iron. Ab-initio data for excitation via 
electron collisions are by far less complete, typically covering transitions up to 
n < 3 or 4 for selected ions of the lighter elements and for iron. Reliable data 
for other members of the iron group and for the heavier elements are only selec- 
tively available. The remainder of data - still the bulk by number - has to be 
approximated for practical applications. 

Consequently, the starting point for the construction of model atoms for many 
elements of astrophysical interest has improved tremendously since the mid-1980s. 
Nevertheless, building realistic model atoms is neither an easy nor a straightfor- 
ward task. It is a common misconception that non-LTE per se brings improve- 
ments over LTE modelling. A careful LTE analysis of well-selected lines can be 
more reliable than a non-LTE study of the 'wrong' lines with an inadequate model 
atom. On the other hand, computations using a realistic model atom will improve 
over LTE - provided that the other ingredients of the modelling are also realistic. 

The independence of microscopic processes from environment - at least under 
not too extreme conditions - provides a tool to assess the quality of model atoms 
by comparison with observation. Comprehensiveness and robustness of a model 
atom are given when it reproduces the observed line spectra over a wide range of 
plasma conditions. Standard stars should serve as test 'laboratories', covering a 
wide range of effective temperature and surface gravity (particle density). 

In the following we discuss practical aspects of the construction of model atoms 
for non-LTE line-formation calculations of trace elements in stellar atmospheres. 
Guidelines and suggestions are given how to build up robust and comprehensive 
model atoms, and how to test them thoroughly. 
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2 Problem Definition 

Non-LTE line-formation calculations solve the coupled statistical equilibrium and 
radiative transfer equations for a prescribed model atmosphere, which may itself 
be in LTE or non-LTE. The computational expenses are therefore lower than for 
full non-LTE calculations that solve also for the atmospheric structure. Hence, 
comprehensive model atoms may be treated in great detail, which turns out to be 
a crucial advantage of restricted non-LTE calculations. 

One of the most important quantities for the comparison with observed spec- 
tral lines - which is the basis of quantitative spectroscopy ~ are accurate occupation 
numbers rii for the levels involved in the transitions. These have to be determined 
in general by solution of the rate equations of statistical equilibrium (though de- 
tailed equilibrium may turn out a valid approximation a posteriori) 

J2 n^{R^J + Qj) = J2 -^J^^n + , (2.1) 

where the Rij and Cij are the radiative and collisional rates, respectively, for the 
transitions from level i to level j. Radiative upward rates are given by 

R^J = At: j n,j{u)^dv , (2.2) 

where crij{v) is an atomic cross-section for bound-bound or bound-free processes, 
the mean intensity, h Planck's constant and v the frequency. In the case of 
collisional processes the upward-rates are given by 



Cij=ne J aij{v)f{v)vdv , (2.3) 

where is the electron density (for the moment we assume that collisions with 
heavy particles can be neglected, see Sect. [5]), v the velocity and f{v) the (in gen- 
eral Maxwellian) velocity distribution of the colliding particles. The corresponding 
downward rates are derived from detailed-balancing arguments, requiring correc- 
tion for stimulated emission in the case of the radiative downward rates. 

Inspection of Eqns. I2.1H2.3I shows which input quantities need to be reliably 
known in order to facilitate an accurate determination of the level populations: 
i) the local temperatures and particle densities, and ii) the non-local radiation 
field, and ill) accurate cross-sections (Jij, and iv) all transitions relevant for the 
problem have to be taken into account. Any shortcomings in these will affect the 
final accuracy that can be obtained in the modelling. 

Items i) and ii) depend on the prescribed model atmosphere used for the re- 
stricted non-LTE calculations. This has to give a fair description of the real tem- 
perature gradient and the density stratification in the star's atmosphere under 
investigation. Particular care has to be invested in the stellar parameter deter- 
mination (effective temperature T^s, surface gravity log 5). This has to account 
for a self-consistent treatment of quantities which are often thought to be of sec- 
ondary nature (microturbulence ^, metallicity, helium abundance, etc.), see Nieva 
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& Przybilla (this volume) . Good knowledge of the model atmosphere code is a pre- 
requisite for successful non-LTE line-formation computations for individual cases 
- use of published model grids allows only to scratch the surface of the problem. 

Items III) and iv) are related to the model atom. Nowadays, the question is 
often not whether to use ab-initio data for the model atom construction, but which 
of the available datasets to adopt. This is a matter of experience and familiarity 
with atomic physics. The quality of agreement of ab-initio results with available 
experimental data should certainly guide the decision. A first important step is 
to check how well the ab-initio calculations reproduce the observed energies of the 
levels in an atom/ion. A comparison of observed and computed oscillator strengths 
and cross-sections for reactions gives further indications (see for example reviews 
by WiUiams 119991 Kjeldsen [5DDBj). The latter will typically be possible only for 
ground and metastable states, but the opportunity for constraining the accuracy 
of the theoretical data should not be missed. Another criterion is the agreement 
between length and velocity forms of oscillator strengths to verify the internal 
consistency of the ab-initio calculations. Eventually, different model atoms may 
be constructed using the alternatives as input data to decide empirically which 
dataset should be preferred for the practical work. 

Experience is also the key in deciding how extended a model atom should be 
and which transitions should be considered. Nature realises all possibilities, but 
we have to handle the mathematical solution of a set of equations describing the 
physics packed into a (restricted) model. In particular, the numerical solution of 
the set of linear equations (j2.1|) requires to be handled carefully. Precautions need 
to be taken to keep the numerical problem well-conditioned and the algorithms 
stable. In practice this means that some transitions may better be ignored in a 
model atom, or the number of levels may be restricted in order to obtain meaning- 
ful results. Larger model atoms accounting for more transitions are therefore not 
necessarily 'better', even if the individual atomic input data are of high quality. 

Finally, execution times of the model computations are also important for the 
practical work. They should not be excessive, requiring a certain compactness of 
model atoms even for non-LTE line-formation computations on prescribed model 
atmospheres. Hence, the construction of model atoms is effectively a highly com- 
plex optimisation problem. A compromise needs to be found between comprehen- 
siveness, accuracy, stability and efficiency. 

3 Model Atom Structure 

The first step in the construction of a model atom is to decide on the extension 
of the model: which ions, which energy levels should be included, and how? This 
depends on the specific non-LTE problem and may range from a few levels for 
studies of a resonance line to many hundred - including packed 'superlevels' - 
if a reproduction of the complex spectra of e.g. iron group species is aspired. 
In the following we concentrate on the general strategy for the construction of 
comprehensive model atoms, which are able to reproduce practically the entire 
observed spectra of an ion over a wide parameter range. 
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Fig. 1. Schematic plot of the Maxwellian velocity distribution for typical temperatures 
at line-formation depths in A- and O-type stars. 
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Fig. 2. Run of departure coefficients bi in a model of the bright supergiant /3 Ori (B8 la) 
as a function of Ti-oss for H I model atoms of different complexity: for 10, 15, 20, 25, 30, 
40, 50 levels (dotted, dashed-dotted, dash-dot-dot-dotted, long dashed, full, dashed, full 
thick lines). The lower right panel displays the behaviour of all levels in the 50-level 
model atom. Note the Rydberg states asymptotically approaching the Hll limit. See the 
text for details. From Przybilla & Butler (UDDH). 
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The fundamental parameter that determines which ions should be included 
in a model atom is the effective temperature, as this determines the energetics 
of the microscopic processes. The term structure of the ions provides a second 
criterion. Hence, the main ionization stage should be adequately represented, plus 
usually two or three minor ionic species - which may comprise in fact the ones 
of particular interest. E.g., the main ionization stage of carbon in a BOV star 
(with Toff 32 000 K) is C ill in the line- formation region, but features of the 
minor species C ii and C IV are also present in the optical spectra. Consequently, 
a comprehensive model atom should consider C ii-iv, plus the ground level of C V. 
The latter has an enormous energy gap between the ground and the first excited 
level, of about 300 eV, such that from atomic structure considerations excited C v 
levels can be safely ignored in the model atom. A lithium model atom for use with 
solar-type stars can be kept simple - detailed Li i -I- the Li il ground level - for the 
same reason, despite Li ii dominating by far in terms of population. 

Concerning the choice of levels to include in a comprehensive model atom there 
are two objective criteria available, which are tightly coupled. One is the energy 
gap between the energetically highest level of the ion and the continuum, as defined 
by the ground level of the next higher ionization stage. As a general rule the gap 
should be less than kT to ensure an accurate determination of ionization fractionqj. 
In order to keep a model atom robust it is recommended to include energetically 
higher levels than the minimum necessary to cope with a given problem. Non-LTE 
studies of the Mgll A4481 A line may serve as an example: the line is observable 
from F-type to mid 0-type stars. For explaining the behaviour of the line in 
early-type stars (Mihalas 11972^ it may be sufficient to consider levels up to n = 5 
or 6 using the criterion above, as the majority of the colliding particles will be 
able to facilitate ionization (see Fig. [T|). However, the same model atom will not 
be useful for analyses of A-type stars, as only a small fraction of electrons in the 
high- velocity tail of the Maxwell distribution will be able to overcome the reaction 
threshold. In that case, completeness up to n = 8 would be required, etc. 

The second criterion is the convergence of the behaviour of level departure co- 
efficients hi with increasing model complexity (for a given set of transition data). 
An example is shown in Fig. [5J for hydrogen model atoms considering levels up 
to n=10 to 50. Models with a too low a number of levels can show a differ- 
ent behaviour at line-formation depths than the more complex models, resulting 
in inaccurate predictions. An alternative formulation of this criterion is via the 
convergence of the line source function (e.g. Sigut ll996| . 

For many elements fine-structure states of a term may safely be combined 
into one level representing the term. This comprises, in particular, the cases that 
are approximated well by LS'-coupling. Collisions couple the individual sub-levels 
tightly because of their small energy separations, i.e. they are in LTE relative to 
each other. A similar opportunity for simplification of the model atom opens up 
for levels at higher excitation energies. The energy separations of terms decrease 



^Minority species are particularly sensitive to non-LTE effects because any small change of 
the ionization rates can affect their populations by a significant amount. 
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Fig. 3. Example of a comprehensive model atom structure. Grotrian diagrams for 
the singlet and triplet spin systems of Mgl are shown. Displayed are the radiative 
bound-bound transitions treated explicitly in the non-LTE calculations. The ionization 
threshold is indicated by the dashed line. From Przybilla et al. l|20Uip . 

with increasing angular quantum number £ for the same n, and in general with 
increasing n. Eventually, they may be safely grouped into one level with appro- 
priate statistical weight. This helps to keep the number of explicit non-LTE levels 
to be treated for elements up to about calcium below ^^200, even if several ion- 
ization stages are treated simultaneously. An example of a comprehensive model 
atom structure is visualised in Fig. [3] in form of a Grotrian diagram, for Mgi. For 
the heavier elements with complex electron structure like the iron group elements 
it may be worthwhile to consider regrouping a multitude of levels with similar 
properties into 'superlevels' (and the transitions into 'superlines'), a concept first 
introduced by Anderson (|1989|) . 

4 Radiative Transitions 

The non-local character of the radiation field drives the stellar atmospheric plasma 
out of detailed equilibrium: photons can travel large distances before interacting 
with the particles, coupling the thermodynamic state of the plasma at different 
depths in the atmosphere. This affects the excitation and ionization of the mate- 
rial. Radiative transitions obey selection rules. 
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Fig. 4. Comparison of oscillator strengths for Cll: values from Froese Fischer & 
Tachiev (|2UU4|) and Nahar POD^ vs. data from the Opacity Project (Yan et al. [1557)1 . 
From Nieva & Przybilla IpOOSl NP08). 

4.1 Line Transitions 

Changes in the internal energetic state of an atom/ion can occur by the absorp- 
tion/emission of photons, giving rise to spectral lines. The strength of a spectral 
line is basically determined by the number of absorbers (or emitters) and the line 
absorption cross-section, which is given by 

2 

where e is the electron charge, TOe the electron mass and fij the oscillator strength. 

is the line absorption profile (the emission line profile is identical assuming 
complete redistribution), which can be approximated well by a Doppler profile in 
most cases when the coupled radiative transfer and statistical equations are solved 
to determine the level population£0 in the restricted non-LTE approach. 

The accuracy of the oscillation strengths used for the model atom construction 
will be a limiting factor for analyses. Fortunately, high-accuracy data is available in 
many casefH, obtained both from experiments as well as from ab-initio calculations. 

It is worthwhile to cross-check available data sources. Newer data are not 
necessarily better, even if a apparently more advanced method was used for their 
determination. As usual, the devil is in the details. An example is shown in Fig.|4l 
Oscillator strengths for C ii from the OP (i?-matrix method in the close-coupling 



■^Detailed line broadening theories become important for tlie last step of non-LTE line- 
formation computations, when synthetic spectra are calculated. 

''Oscillator strengths are the only ones among the atomic data discussed further on in Sects. |4]- 
[6] that are required in LTE investigations. 
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Fig. 5. Comparison of C ll pliotoionization cross-sections from tlie Opacity Project (Yan 
& Seaton IT987|) and Nahar & Pradhan ([T997)l for and 4/ as a function of 

wavelengtli. From NP08. 



approximation assuming L5-coupling) are compared to data of Nahar (|2002p ob- 
tained with the (in principle superior) Breit-Pauli _R-matrix (BPRM) method and 
the smaller number of high-precision data from application of the Multiconfigura- 
tion Hartree-Fock (MCHF) method. The in general good agreement between the 
OP and the MCHF data (with few exceptions they follow the 1:1 relation) and the 
considerable scatter in the comparison of OP and Nahar's /-values was one of the 
factors to disregard this particular set of BPRM data from the model atom con- 
struction in that case. However, the question which of the available data should 
be used for the construction of model atoms has no simple answer in general. 

Limitations on the number of transitions that can be considered explicitly in the 
radiative transfer calculations may be imposed by the numerical method chosen 
for solving the equation systems. Complete linearisation techniques are much 
more restricted than Accelerated Lambda Iteration (ALI) techniques. The most 
important transitions (high /-value, location at wavelengths with non-negligible 
flux) will need to be considered in the former case while virtually all line transitions 
may be accounted for in the latter case, see Fig. [3] for an example. 

4.2 Photoionizations 

The availability of photoionization cross-sections has improved enormously in the 
past 20 years, mostly because of the efforts made in the OP and IP. Like for oscil- 
lator strengths, the quality of the different data sets should be evaluated for the 
model atom construction. A comparison of photoionization cross-sections from 
two different calculations is shown in Fig. [SJ for two excited levels in C ii. Maybe 
surprisingly, some details can be relevant for the model atom construction, while 
obvious discrepancies may be not. The small shifts in the resonance structure 
near the ionization threshold of the energetically relatively low- lying 2p^ state 
of C II (see Fig. [8] for a Grotrian diagram) can be important for the rate determi- 
nation, as they occur near the flux maximum of OB-type stars. The more than 
two orders-of-magnitude differences in the high-energy tail of the cross-section for 
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the 4/^F° level (occurring because of the different targets in the two ab-initio 
calculations) are irrelevant on the other hand, because of the low flux at these 
short wavelengths. It has to be decided by comparison of the results from differ- 
ent model atom realisations which of the available data should finally be selected, 
see Sect. [7] for an example. As always, some guidance may also be obtained by 
comparison of the different theories with experiments. 

The OP and IP have provided ab-initio data for photoionization from levels 
with typically n < 10 and £ ~ 4. Missing data for levels of higher n or £ may be 
rather safely approximated by hydrogenic cross-sections (Mihalas ll978| p. 99) 

aH= 2.815x10^9^^11%^, (4.2) 

for the model atom construction, in particular as these levels are usually packed, 
see Sect. [3l Here, Z denotes the charge of the ion and gn the bound-free Gaunt 
factor, which is of order unity at ionization threshold. The threshold cross-section 
is given by (^(i^o,^^) = 7.91 x 10^^^ /Z^ ngu cm^. 



5 Collisional Transitions 



Inelastic collisions with particles can also lead to excitation and ionization of 
atoms/ions. The velocity distribution of particles in stellar atmospheres is in prac- 
tically all cases Maxwellian, determined by the local plasma temperature. Collision 
processes will therefore drive the plasma towards LTE. Contrary to radiative tran- 
sitions, no selection rules apply. Typically, only electron collisions are considered, 
as the thermal velocity and hence the collision frequency of heavy particles is much 
smaller, e.g. by a factor ^-^43 for hydrogen. This is in particular valid for all envi- 
ronments where the stellar plasma is sufhciently ionized. Heavy particle collisions 
may become important in special cases, however, like for cool metal-poor stars. 



5.1 Collisional Excitation 

The main aspect of collisional excitation by electron impact is certainly the ther- 
malising effect on level populations in general, dampening out departures from 
detailed equilibrium imposed by the radiative processes. Curiously enough, col- 
lisional excitation can also drive individual level occupations out of LTE under 
special circumstances. This occurs for cases where a level is collisionally tightly- 
coupled to another level that experiences strong non-LTE departures. Examples 
are energetically close levels from different spin systems, one being metastable (a 
radiative decay to the ground state is prohibited by selection rules) . A prominent 
example for this coupling are the 3s '^S and 3s ^S levels of O i (Przybilla et al. 120001 
Fabbian et al lMig]) . 

Larger sets of collisional excitation data for transitions up to typically n = 3 or 
4 are available now, e.g. from the IP. A good part of the data are published in 
the astrophysics literature, but the reader should be aware that much more data 
can be found in the physics literature because of their relevance for fusion research 
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Fig. 6. Examples of collision strengths Q for optically allowed and forbidden transitions 
in C II in the near-threshold region as a function of incident electron energy (fine-structure 
data from Wilson et al. I2005p . For comparison, = 1 is indicated by a dotted line. 



and technical applications. Most useful for practical applications are tabulations 
of thermally-averaged effective collision strengths 



POO 

Tij^ / nije^p{-Ej/kT)d{Ej/kT) , 
Jo 



(5.1) 



where Ej is the energy of the outgoing electron and fiij the collision strength, which 
is symmetric as well as dimensionless, ftij = flji. Effective collision strengths 
facilitate an easy evaluation of transition rates, which are proportional to Tij. 

While ab-initio data are of highest relevance for the construction of model 
atoms, one has to resort to approximation formulae for the vast majority of possible 
transitions. Different descriptions are available from the literature, a comparison 
can be found in Mashonkina (|1996p . In the following we concentrate on the two 
most commonly used approximations. 

The semiempirical formula of Van Regemorter ()1962p allows rates for radia- 
tively permitted transitions to be evaluated in terms of the oscillator strength. 



a, = 5.465 X 10-"neri/2[14.5/„(4/uo)]exp(-uo)r(uo) : 



(5.2) 



where mh is the ionization potential of hydrogen scaled by kT, uq = E^/kT [Eq 
is the threshold energy for the process), and r(uo) = max[g, 0.276 exp(uo)i?i (uo)] 
for ions (Ei is the first exponential integral). The parameter g is about 0.2 if the 
principal quantum number changes during the transition and about 0.7 if not. For 
neutral atoms T{uo) has a different form (see Auer & Mihalas [1973]) . 
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Excitation rates of radiatively forbidden transitions are often evaluated accord- 
ing to the formula of Allen (|1973p . 



where Fy (T) is a temperature-dependent factor in the collision rate as defined by 
Mihalas (|1978l p. 133). Typically, 51 = 1 is assumed for forbidden transitions. 

Note that Eqns. [5^ and [5731 give, at best, order-of-magnitude estimates. Com- 
parisons of ab-initio data with the approximations should be made whenever possi- 
ble to evaluate the true uncertainties that can be expected. An example is shown 
in Fig. ini Differences of up to several orders of magnitude are found, and the 
quantum-mechanical data can show pronounced resonance structure. It is there- 
fore advisable to investigate the available ab-initio data for trends and regularities. 
Extrapolations can be made based on these in order to improve the approximate 
input data for the model atom construction. 



5.2 Collisional Ionization 

The rates of collisional ionization by electron impact are mostly affected by the 
availability of electrons with energies high enough to overcome the threshold for the 
reaction. This makes collisional ionization from the ground state and energetically 
low-lying levels rather inefficient, as only few electrons in the high-velocity tail of 
the Maxwell distribution are available for this. On the other hand, collisions 
become a dominant factor for the coupling of high-lying levels to the continuum. 

Unfortunately, cross-sections for ionization by electron impact are among the 
least-constrained atomic data. Experiments usually cover ionization from the 
ground state only, and the reader is referred to the atomic physics literature for 
extracting data for a specific problem. On the theoretical side, few methods have 
been successfully applied. The fundamental challenge which distinguishes colli- 
sional ionization from excitation is the fact that the Coulomb interaction between 
each of the two outgoing electrons and the residual ion is present even at large dis- 
tances. Recently, breakthrough results have been obtained by use of the converging 
close-coupling method. Several fundamental processes have been modelled accu- 
rately, providing cross sections that closely reproduce the available experimental 
data in these cases, see Bray et al. (|2002p for a review. 

However, in the majority of cases one has to rely on more simple approaches. 
An often used approximation formula for quantifying collisional ionization was 
given by Seaton (|1962p . which expresses the collisional cross-section in terms of 
the photoionization cross-section, yielding a rate 

= 1.55 X 10^^neT^^/^g,a{i^o) exp(-iio)/uo , (5.4) 

where (t(i'o) is the threshold photoionization cross-section (Sect. 14. 2p . and gi is 
of the order 0.1, 0.2, and 0.3 for Z = l, 2, and >3, respectively. Again, Eqn. 15.41 
provides, at best, an order-of-magnitude estimate. A comparison of collision rates 
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calculated with Eqn. l5.4l witli those evaluated by an empirical analytical expression 
(Drawin I196ip indicates that the uncertainties may be sometimes much larger 
(^Mashonkina 11995)) . The Seaton formula (and analogous simple approximations) 
should therefore be applied with caution. 

5.3 Hydrogen Collisions 

Excitation and ionization by inelastic heavy particle collisions are usually consid- 
ered unimportant in comparison to electron collisions, which occur more frequently. 
However, the ratio of hydrogen atoms to electrons may easily exceed a factor 10^ 
in cool stars, in particular in metal-poor objects. In such a case H collisions may 
become a dominant thermalisation process, which must not be neglected. 

Characteristic collision energies in cool star photospheres are of the order 
kT w 0.2 — 0.6 eV, i.e. they are comparable to typical atomic transition ener- 
gies. Consequently, the near-threshold behaviour of the cross-sections is most 
important for the determination of the collision rates. Laboratory measurements 
or ab-initio calculations of cross-sections near threshold are scarcely found in lit- 
erature for these neutral particle collisions. Some progress has been made recently 
for the Na + H system (Fleck et al. [T99T| Belyaev et a. [T9991 [20T0| . showing that 
the collision rates can differ by several orders of magnitude compared to simple 
approximations (see the discussion below). Similar discrepancies were also found 
for the Li -I- H system (Belyaev k Barklem [2003|) . 

In face of the absence of reliable data for practically all cases of interest, one 
has to resort to the use of approximations for the description of inelastic hydrogen 
collisions. Most work relies on the Steenbock & Holweger (|1984p approximation, 
who generalised Drawin's formulae (Drawin ll9681fT969l based on Thomson's clas- 
sical theory), originally developed for collisions between identical particles. The 
Maxwell-averaged formulae for excitation and ionization of particle species A by 
collision with H becomes 



with /i — toatoh/(wa + ™h) denoting the reduced mass. For coUisional exci- 
tation Q — AE)"^ flu and W = AE/kT and for collisional ionization Q = 
{In/ Ia)^ fi£,i and W = /a/A:T, where AE denotes the energy difference between 
upper and lower state of the transition and /h and /a are the ionization energy 
of hydrogen and atomic species A; fiu is the oscillator strength, fi is an efficient 

oscillator strength for ionization and the number of equivalent electrons. The 

2 2m 
function ^'(M/) is given by *(H/) = (l + — )exp(-W)(l + - ^tttt)"^- Note 

that Eqn. l5.5l does not apply to optically forbidden transitions. Takeda (I199ip sug- 
gested to relate the hydrogen collision rate with the electron collision rate via 




(5.5) 




(5.6) 
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assuming ad-hoc a similarity of cross-sections for both cases {Cfj is usually evalu- 
ated with f2 = 1). Here, njj denotes the number density of neutral hydrogen. 

Often the results are scaled by a factor ^h- A value of S'h = is equivalent 
to no hydrogen collisions, S'h oo enforces LTE. In general, S'h is constrained 
empirically by demanding that the scatter in abundance as determined from the 
entire sample of lines of a species should be minimised. 

Equations 15.51 and 15.61 were considered appropriate to provide an order-of- 
magnitude estimate for a long time, but in view of the few ab-initio calculations 
available now, this assessment appears too optimistic. This confirms earlier indi- 
cations of an underestimation of the real uncertainties, from the comparison with 
other approximations, see e.g. Mashonkina (|1996p . Barklem (|2007p investigated 
the uncertainties for the most simple case, H -|- H, in detail. 

The Steenbock & Holweger formula remains in use for astrophysical applica- 
tions in view of the lack of other reliable data, despite all the warning evidence. 
In view of this, efforts should be made to determine proper scaling factors Sh in 
order to minimise the impact on the accuracy of analyses. This could be made by 
extensive investigations following the recommendations for testing model atoms 
in Sect. [71 covering the parameter space comprehensively (wide range of effective 
temperature, densities and metallicity) . 

6 Other Processes 

Consideration of the processes described in the previous two sections is usually 
sufficient for the construction of model atoms for the analysis of stellar photo- 
spheres. Nonetheless, two other types of radiative and collisional processes may 
be of relevance in some cases. They are only briefly described in the following for 
completeness, leaving it to the reader to investigate the specialist literature. 

Autoionization can occur in complex atoms with several electrons. When two 
electrons are excited simultaneously, this can give rise to states below and above 
the ionization potential. The states above the ionization threshold may autoionize 
to the ground state of the ion, releasing one electron. The inverse process is also 
possible, and, if a stabilising radiative decay occurs within the (short) lifetime of 
the doubly excited stat^, it can give rise to an efficient recombination mechanism. 
This is referred to as dielectronic recombination. Details of rate coefficient mod- 
elling can be found e.g. in Badnell et al. (|2003p . An application in the context of 
WR-type central stars of planetary nebulae is discussed by de Marco et al. (|1998p . 

Charge exchange reactions are collisional processes between atoms/ions in which 
one, or more, electrons are transferred from one collision partner to the other, 
e.g. A-t-B"*" ^> A+ -f B, with B usually being H or He. One well-known reaction 
occurring in stellar atmospheres is O -I- H+ ^ 0+ + H, which can dominate the 
ionization balance of oxygen as the non-LTE departures of the 7i(Hi)/n(Hii) are 



*The presence of these so-called autoionizing states has consequences for the absorption of 
photons and the scattering of electrons by atoms/ions: it gives rise to resonances in the pho- 
toionization and electron collision cross-sections, see Figs. [5]and|6l 
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Fig. 7. Left: sensitivity of two Cll lines to use of photoionization cross-sections from 
different ab-initio calculations in a model spectrum for a B2 V star. Right: differences in 
(logarithmic) carbon abundance as derived from the two lines in 6 stars as a function of 
effective temperature. From NP08. 



forced upon n(0i)/n(0ii) by this resonant reaction (the ionization potentials of 
neutral hydrogen and oxygen are very similar), see e.g. Baschek et al. ()1977p or 
Przybilla et al. (|2000p for a discussion. A list of reaction playing a possible role in 
astrophysical plasmas and tabulated reaction rates are provided e.g. by Arnaud 

6 Rothenflug (OMSll . 

7 Testing Model Atoms 

The question whether a model atom is realistic can only be answered by compari- 
son with observation. One needs to test whether the model atom is comprehensive 
enough, i.e. whether the level structure and all relevant connecting transitions are 
considered properly and whether the atomic data used are sufficiently accurate. 
Usually, this should give rise to an iterative process: a stepwise improvement of the 
model atom by empirical selection of the 'best' input data. The aim is to single out 
one model atom, that reproduces the observed spectra closely at once, independent 
of the plasma conditions (atomic properties are independent of environment). 

In order to give an idea on the practical approach for performing such tests 
we discuss an example. Synthetic line profiles from calculations with two different 
model atoms are compared in the left panel of Fig. [71 While the strong C ii 
A4267A transition is highly sensitive to non-LTE effects - in particular to the 
photoionization cross-sections adopted the other line is virtually insensitive to 
any model atom realisation using reasonable atomic input data (it is 'in LTE'). 
Such a sensitivity is one of the keys to select the 'right' photoionization data for the 
model atom construction. The second ingredient in this process is the comparison 
with observations, here for stars in a temperature sequence (right panel of Fig. [7])- 
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Fig. 8. Grotrian diagram for the Cll doublet and quartet spin systems. Transitions that 
give rise to lines in the optical spectra of B-type stars are identified by continuous and 
dotted lines. These are the observables available for testing model atoms. Additional 
channels become available in the UV and IR spectral range. From NP08. 

This is in order to test the reaction of the model realisations to a hardening 
radiation field. The line 'in LTE' serves as reference and the goal of the model 
atom optimisation is to minimise the differences in the derived abundances from 
the various indicators. In this case it was shown that ill-chosen atomic data can 
result in line abundance differing by up to 0.8 dex (NP08), which helped to solve 
one notorious non-LTE problem of stellar astrophysics (Nieva & Przvbilla l2006p . 

Ideally, the observations used for the model atom calibration should range from 
the far-UV to the near-IR in order to test as many channels (spectral lines) as possi- 
ble, even those that may be irrelevant for practical applications later. At the same 
time, the observations should cover an as wide range of plasma parameters as possi- 
ble: high-density, i.e. collisionally dominated, environments like the photospheres 
of dwarf stars and low-density plasmas (those dominated by radiative processes) 
as encountered in (super-)giants should be considered, at different temperatures. 
Where possible, also the metallicity dependence of non-LTE effects should be in- 
vestigated to test the response of a model atom to a varying radiation field, e.g. 
by considering stars of Population I and II. Further tests may include more 'ex- 
otic' environments, like He-dominated plasmas in compact subdwarfs (Przybilla 
et al. I2006a|) or giant extreme helium stars (Przybilla et al. .2005, .2006b)) . 
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Such a comprehensive approach involving satellite observations may not be fea- 
sible in almost all cases. However, high-quality spectroscopy from the ground with 
modern echelle spectrographs is often sufficient to provide the means to facilitate 
proper model atom testing. Figure [5] visualises the channels available for testing 
a model atom of C II using the optical spectra of early B-type stars. Note that 
despite this comprises a fair number of energy levels and transitions from several 
multiplets, only a small fraction of the entire model atom can be really scrutinised 
by this. Consideration of lines from additional ionization stages (C ill and C IV 
in that case, NP08) may put further constraints, as the full set of observed lines 
should be reproduced closely by the model simultaneously. 

The above example is typical for elements with relatively simple electron struc- 
ture. More complex electron configurations like the open Sd-shells in the iron-group 
members pose a larger challenge at first glance, but the enormous number of ob- 
servable transitions puts many constraints on the model atom construction as well. 
The real challenge are therefore the low-abundance light elements lithium, beryl- 
lium and boron, and in particular their alkali-like ionqj. There, typically only the 
resonance lines are observable, which gives only marginal constraints for tests of 
the model atom. Where possible, resonance lines from another ionization stage 
should therefore be investigated, or subordinate lines that may become observable 
in stars with particular high abundance in the element under study. 

It is of utmost importance to use realistic atmospheric structures for testing 
model atoms, requiring a proper determination of the stellar atmospheric parame- 
ters, see Nieva & Przybilla (this volume). Well-studied standard stars like the Sun 
(G2V), Procyon (F5IV-V), Vega (AOV), rSco (B0.2V) or Arcturus (K1.5III) 
with tight constraints on their atmospheric parameters are therefore primary ob- 
jects for the comparison of the models with observation. However, further stars 
that bracket the extremes of the parameter space to be studied have also to be 
considered. 

8 Final Remarks 

The plethora of possibilities implies that the model atom construction does not 
result in a unique solution. A good reproduction of observations may be achieved 
by a whole family of models. The main insight is that there are many insufficient 
model atoms, but few adequate ones. In consequence, non-LTE analyses are not 
superior to LTE investigations per se, but require robust and comprehensive model 
atoms. In view of the increasing availability of accurate and precise atomic data 
from ab-initio calculations and experiments one is faced by a perpetual challenge: 
the impact of new high-quality atomic data should be tested on the modelling 
whenever such become available. Of course, the same is true for new observations 
that may facilitate the predictive power of the models to be tested further by 
opening up other channels for the model atom calibration. 



^High-quality atomic data can fortunately be obtained with relative ease for these simple ions. 
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